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Abstract The angular distribution of electrons accelerated in solar flares is a 
key parameter in the understanding of the acceleration and propagation mech- 
anisms that occur there. However, the anisotropy of energetic electrons is still 
a poorly known quantity, with observational studies producing evidence for an 
isotropic distribution and theoretical models mainly considering the strongly 
beamed case. We use the effect of photosphcric albedo to infer the pitch angle 
distribution of X-ray emitting electrons using Hard X-ray data from RHESSI. A 
bi-directional approximation is applied and a regularized inversion is performed 
for eight large flare events to deduce the electron spectra in both downward 
(towards the photosphere) and upward (away from the photosphere) directions. 
The electron spectra and the electron anisotropy ratios are calculated for broad 
43 ■ energy range from about 10 and up to ~ 300 kcV near the peak of the flares. 

The variation of electron anisotropy over short periods of time intervals lasting 
4, 8 and 16 seconds near the impulsive peak has been examined. The results 
show little evidence for strong anisotropy and the mean electron flux spectra are 
consistent with the isotropic electron distribution. The 3<7-level uncertainties, 
although energy and event dependent, are found to suggest that anisotropic 
distribution with anisotropy larger than ~ 3 are not consistent with the hard 
X-ray data. At energies above 150 — 200 kcV, the uncertainties are larger and 
thus the possible electron anisotropics could be larger. 

Keywords: Flares; X-Ray Bursts, Spectrum; Energetic Particles, Electrons; 
Corona 



1. Introduction 

Solar flares are one of the most energetic processes which occur in the solar sys- 
tem. However the details of the acceleration processes responsible are still poorly 
known. The X-rays which are observed at the Earth are commonly produced 
by the accelerated electrons interacting with the solar plasma and producing 
bremsstrahlung radiation. To understand solar flares it is therefore important 
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to understand the distribution of these accelerated electrons. In general, this 
distribution will vary in space, energy and in pitch angle (see IHohnan et al.i 
120111 IKontar et al.i 120111 as the recent reviews of electron properties in solar 
flares). 

Several techniques have been used to estimate the anisotropy in the pitch- 
angle distribution of X-ray emitting electrons in solar flares (Holman et ft/.. 2011; 



Kontar et al., 2011). The most commonly used method is to look at the centre- 



to-limb variation of solar X-ray properties ( |Datlowe, Elcan, and Hudson, 1974 
|Vestrand et al, 1987fl . That is comparing the characteristics, most commonly to- 
tal X-ray flux (Pizzichini, Spizzichino, and Vcspigna ni, 1974| , or the spectral in- 
dex of solar flares at the limb to disk centre events. Studies concentrating on lower 
energy emission (below 300 keV) tended to find no significant evidence of direc- 
tivity. These studies have also been performed using SMM data ( |Dennis, 1988 ) 
studying flares with energies above 300 keV ( |Vilmer,T9 94; Bai, 1988; Bogova lov et al., 1985p .[ 
Some evidence for directivity at high energies has been reported ( McTicrnan an d Petrosian, 1991| [ 



Vestrand, Forrest, and Rieger, 1991 1. More recently RHESSI (Ramaty High En- 



ergy Solar Spectroscopic Imager) data ( |Lin et al, 2002[ ) has been used to deter- 
mine the X-ray anisotropy ( |Kasparova, Kontar, and Brown, 2007| ). An obvious 
disadvantage of the statistical method is that the variation can only be seen as 
an average over a large number of solar flares, so little can be said about X-ray 
or electron anisotropy in a given flare. 

An approach which allows individual flares to be studied is the stereoscopic 
method (jCatalaii o and van Allen, 1973[ ). Here each individual flare is measured 
directly by two spacecraft at two different locations, ideally well separated in 
space. Studies that have been performed using this method do not show any 
clear evidence of directivity ( |Kane et al, 1998] |Li et al, 1994[ ) . A disadvantage 
of this approach is the difficulty in cross calibrating, often leading to large errors. 
Another drawback is, as in the centre-to-limb method, this technique does not 
give direct information about the downward electron distribution. 

As an anisotropic electron distribution will produce polarised X-rays, mea- 
suring the polarisation can therefore give a measure of the anisotropy of the 



electron distribution ( Leach and Petrosian, 1983 ). An isotropic source should 

show low polarisation, whereas a beam should produce significant polarisation 

( |Bai and Ramaty, 19"78l|Leach and Petrosian, 1983||Emslie, Bradsher, and McConnell, 2008[ ).| 

Studies of polarisation have been performed using various X-ray satellites, re- 



cently using the Coronas-F satellite ( Zhitnik et al., 2006 ). The reported mea- 
surements vary substantially from observation to observation, adding to the 
scepticism of these measurements. A major drawback of this approach is the 
observational difficulty in measuring polarisation at HXR energies for tran- 
sient events like solar flares. Several attempts have been made using RHESSI 
(Mc Connell et al., 2002 Suarez-Garcia et al., 2006 ), but so far there have been 
no conclusive measurements made. HXR polarisation has not yet been used to 
its full potential and future observations could provide a more definitive answer. 

Another important process, which can be used to diagnose the angular distri- 
bution, is photospheric albedo ( |Kontar and Brown, 2 006). Photospheric albedo 
of X-rays results from initially downward directed X-rays, which are Comp- 
ton scattered off electrons in the solar photosphere, being observed at Earth. 
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This effect distorts the HXR spectrum QTomblin, 1972] |Bai and Ramaty, 1978| 
and leads to an appearance of apparent low energy cutoff in the deduced elec- 
tron spectrum ( jKontar, Dickson, and Kasparova, 2008 ), when the observed X- 
ray spectrum is not corrected for an albedo ( Kontar et ai, 2006] ) . As the spectral 
shapes of reflected and primary hard X-ray spectra are sufficiently distinct, these 
two components can be distinguished and the albedo reflected flux could be 
used as a measure of the downward going electrons (Figure 1). RHESSI provides 
sufficient energy resolution, broad energy coverage, and sensitivity to better 
constrain directivity of energetic electrons in individual solar flare events. 

Here wc use this albedo method to examine the directivity of energetic elec- 
trons in solar flares. We first perform some forward modelling to examine the 
effect directivity of the electron spectrum and the role of the albedo contribu- 
tion has on solar flare spectra. The RHESSI flare catalogue has been searched 
for suitable flares between 2002 and 2008 and we use the spectral data from 
the impulsive phases of several well observed flares to perform a bi-directional 
inversion, estimating the fluxes of electrons travelling towards and away from 
the photosphere. 



2. Solar flare spectrum and electron anisotropy 

The spatially integrated X-ray photon flux spectrum at the Earth /(e, 9q), [phot- 
ons cm -2 s _1 kcV" 1 ] is straightforwardly related to the mean electron flux 
spectrum of energetic electron F(E,f3) [electrons cm~ 2 s _1 kcV -1 sr^ 1 ] via the 
linear integral relation (e.g. IMassone et al.\ (|2004|) . |Jeffrey and Kontar] (|201ip ). 
For a given electron distribution F(E,ft), where ft is the electron pitch angle, 
the resulting emitted X-ray flux at distance R is 

/(£ ' 6o) = iS* fj f T{E > /?)Q(e ' E > ^)) sm (P)dEdftd^ , (1) 

where V is the source volume, h is the mean plasma density, Q(e,E,9) is the 
angular dependant bremsstrahlung cross-section, e is photon energy, E is electron 
energy and 8 is the angle between the initial electron velocity vector and the 
direction of the emitted photon. The cross-section used here is the electron-ion 
cross section, formula 2BN from lKoch an d Motz ( 1959) with Coulomb correction 
bv lElwertl (|1939[) added. The relation between the angles ft, f?o, (f> and 8 is given 

by 

cosd = cos ft cos8q + sinft sinOo cos<f> , (2) 

where 6$ is the angle between the emitted photon and the direction of the photo- 
sphere. The exact plasma density distribution n and flaring volume V are often 
poorly known and therefore it is preferable if the value nVF(E, fi) is inferred, 
this is the mean electron flux spectrum ( |Brown, Emslie, and Kontar, 2003 ), a 



density weighted electron flux, multiplied by nV the number of electrons in the 
emitting volume. This value is model independent, and its detailed structure can 
be related to the electron acceleration and propagation physics which is central 
to the understanding of solar flares. 
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Figure 1. The geometry of the X-ray emitting source above the photosphere and bi-direc- 
tional approximation. X-rays are emitted in all directions and observed directly at Earth 
or Compton back-scattered in the solar photosphere and then observed at Earth. The true 
angular distribution of electrons F(E, /3) is approximated by downward and upward F u 
going electrons. 

The electron distribution of emitting electrons is assumed to be separable in 
energy and angle. The energy dependence is taken to be a single power-law as 
solar flare x-ray observations are often well fit by power laws implying close to 
power-law electron spectra. For simplicity the pitch angle distribution is taken to 
be a gaussian beam centred downwards. This pitch angle distribution is assumed 
as electrons tend to stream along magnetic field lines resulting in a downwards 
beam, and the gaussian parametrization of this is popular as it has useful analytic 
properties. The total distribution has the form similar to ILeach and Petrosianl 



where /i = cos([3) and Afj, characterises the width of the Gaussian. It should 
be noted that Equation (3) refers to the mean electron flux spectrum not the 
injected electron spectrum considered in ILeach and P etrosianl [19831 Three dif- 
ferent cases were assumed for the angular variation of the electron spectrum 
(Figure 2): the isotropic case (Afj, = 10), an intermediate anisotropic case 
(Afj, = 0.4) and a highly beamed case (Afj = 0.1). The electron spectral index 
S was assumed to be S = 2 as the typical flare spectral index of X-rays in large 
flares is around 7 = 3 (see e.g. IDennisl I1985P and for the mean electron flux 
spectrum 7 » 5 + 1. 

As the assumed electron distribution has azimuthal symmetry an average 
cross-section integrated over <f> can immediately be defined leaving only the an- 
gles 6*o and (3 needed to characterise the angular distribution (c.f. IMassone et al\ 



J 0=0 

The angular distribution of the emitted X-rays with respect to the downward 
direction can now be found by applying this cross section to the assumed electron 



(IT98311 




(3) 




(4) 
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spectrum and integrating over electron energy and pitch angle (Figure 3) 

nv r f°° — 

J ( £ '^) = 4^2 y o J F(E,ft)Q'(e,E,ft,6 )sin^)dEdft . (5) 

The primary emission observed for a flare at heliocentric angle 8' is expected 
to come from a small range in angle in the direction of the observer. Considering 
the geometry, it is clear that the upward directed photon distribution can be 
approximated as 

I u (e)=I(e > = 180 o -9') . (6) 

The albedo reflected component, on the other hand, results from the photons 
directed down towards the photosphere. This is likely to be a broader distribution 
so an average is taken over a downwards directed cone concentric with the mean 
direction of the electron distribution and with half angle a, in this study a value 
of 90° is used. The downward directed flux can then be defined as 

I(e, 0o ) sm 0o^o 

sin0o<i0o 



Id(c) = J -^^ . (7) 



i=0 

The reflected component, due to albedo, of a given X-ray spectrum incident on 
the photosphere can be characterised by using a Green's function A dependant 
on observation angle 0' ( |Kontar and Brown, 2006] ) . Thus, the total observed X- 
ray spectrum will be given by the sum of the directly observed and reflected 
components: 

Io = Iu + AI D , (8) 

where Io is the total and Ijj, Id are upward and downward directed components 
(Figure 6). 

In practice the photon spectrum is calculated numerically with a finite resolu- 
tion. For energy a pseudo-logarithmic binning scheme with 100 bins starting at 
10 kcV and going up to 5 MeV is used for both electrons and photons. Due to the 
highly anisotropic nature of the cross-section at high energies fine resolution in 
angle was needed. Angles ft and 0o were binned in 90 evenly spaced bins between 
and 7r whereas </> was binned in 180 evenly spaced bins between and 2tt. 

The upward and downward components of the photon flux are now given by 
vectors. The observer directed component is defined as 

Iu = {Iu(ei), —, Iu(e n )}, i = l,..,n (9) 

where tj corresponds to the centre energy of the photon bin and n is the number 
of bins in photons space. The photosphere directed component is given by Id 
defined in the same way. Similarly the angle dependant Green's function is here 
calculated in the form of a n x n Green's matrix A. The vector representing the 
total observed flux is simply given by 

Io = Iu + AI D , (10) 
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where A is the albedo matrix constructed from Green's functions. 
2.1. Photon Spectral Index 

Due to the effect of albedo one of the most notable variations with changing 
anisotropy is in the photon spectral index, the derivative with respect to en- 



ergy of the photon flux. It is commonly defined as (Brown and Emslie, 1988 
|Conway et al, 2003D 

. . e dl dlogl 
I de d log e 

which is calculated for both the primary photon spectrum for a range of view- 
ing angles (Figure 5) and for the total observed spectrum including albedo 
component for a range of positions on the solar disk (Figure 7). 




o.o 

1.0 -0.5 



Figure 2. Polar diagrams of the assumed electron pitch-angle distribution F(E, /3). The angle 
made with the x-axis corresponds to the pitch angle /? of the electron and the radial extent 
corresponds to the magnitude of the distribution. The electron flux with fi, = 1 is normalised 
to 1. Left: isotropic case Afi = 10. Middle: intermediate anisotropic case A^t = 0.4. Right: 
Highly beamed case Afi = 0.1. 




Figure 3. Polar diagrams of the emitted photon distribution I(e,8o). The angle made with 
the x-axis corresponds to the angle 8o and the radial extent corresponds to the magnitude of 
the photon distribution. The energy distribution is plotted for several energies: 10 keV (solid 
black), 40 keV (dotted purple), 150 keV (dashed green), 600 keV (dashed yellow) and 5 MeV 
(solid red). The photon flux with cos$o = 1 is normalised to 1, Left: isotropic case Middle: 
intermediate anisotropic case. Right: Highly beamed case. 



This method was performed for the cases of strong (A/z = 0.1) and intermedi- 
ate (A/z = 0.4) beaming over an energy range 10 keV to 5 MeV (Figure 2). After 
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Energy (keV) Energy (keV) Energy (keV) 

Figure 4. Photon flux for several selected values of 0n - 0° (red, long dashed), 45° (yellow, 
dot-dashed), 90° (green, dashed), 135° (blue, dotted) and 180° (black, solid). Left: isotropic 
case Middle: intermediate anisotropic case. Right: Highly beamed case. 




100 1000 
Energy (keV) 

Figure 5. Photon spectral index for several selected values of 8q - 0° (red, long dashed), 45° 
(yellow, dot-dashed), 90° (green, dashed), 135° (blue, dotted) and 180° (black, solid). Left: 
isotropic case Middle: intermediate anisotropic case. Right: Highly beamed case. 



applying this assumed electron spectrum to the bremsstrahlung cross-section 
the angular dependant X-ray emission is found (Figure 3), for lower energies 
this tends towards being closer to isotropic than the electron distribution but 
for high energies it is reasonably similar to the input electron spectrum. The 
emitted flux density (Figure 4) and spectral index (Figure 5) for a range of angles 
of observation are then calculated. A Green's matrix corresponding to the albedo 
reflection at a range of heliocentric angles ( Kontar et al, 2006 1 was then applied 
(Figure 6). Emission close to the solar limb shows very little influence from albedo 
as expected and the power law in photon energy is recovered, however emission 
closer to the disk centre shows a distinctive hump over the entire energy range 
due to the albedo reflection. The influence of albedo can be seen more clearly 
when observed photon spectral index is considered (Figure 7) , flares close to the 
disk centre show a large increase in 7 above 200 keV, and this is more pronounced 
in the strong beaming case. 



3. Application to RHESSI data 

The RHESSI data archive was examined for flares with emission above 300 keV 
with particular attention paid to flares close to the solar disk centre. These flares 
are selected because the forward modelling suggests that the variation due to 
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Energy CkeV) Energy CkeV) Energy (keV) 

Figure 6. Total observed photon flux including reflected albedo component for flares located 
at different places on the disk ranging from disk centre cos 9' = 1 (black), cos 9' = 0.9 (purple), 
cos 6' = 0.8 (indigo) cos6»' = 1 to limb cos0' = 0.2 (yellow) cosff' = 0.1 (orange) cos6»' = 0.01 
(red). Left: isotropic case Middle: intermediate anisotropic case. Right: Highly beamed case. 




lO lOO 10 loo io loo 

Energy (keV) Energy (keV) Energy (keV) 

Figure 7. Photon spectral index for the total observed spectrum against photon energy for 
flares located at different places on the disk ranging from disk centre cos 9' = 1 (black) , 
cosO' = 0.9 (purple), cos6»' = 0.8 (indigo) cos0' = 1 to limb cos6»' = 0.2 (yellow) cos 9' = 0.1 
(orange) cos 9' = 0.01 (red). Left: isotropic case Middle: intermediate anisotropic case. Right: 
Highly beamed case. 



beaming is strongest at high energies and flares closest to the disk centre should 
have the strongest albedo reflection and should therefore also show the greatest 
change due to beaming. In total eight suitable flares were found (Table 1) which 
were within 60° of the solar centre (Figure 8) and showed significant > 3er counts 
above background (Figure 9); a number of other flares matching these criteria 
were found but they were discounted due to high levels of pulse pilcup and 
particle contamination. 

For each of the flares found the background was removed in the standard 
manner (Schwartz et al., 2002). Counts were accumulated over the impulsive 
phase, as the differences in the spectra due to anisotropy in the electron spectra 
are greater at higher energies. The time intervals studied were selected ensuring 
a high number of high energy counts (Figure 10). A pseudo-logarithmic binning 
scheme between 10 keV and 500 keV was used to initially accumulate the spectra 
avoiding detectors 2 and 7 due to their poor resolution ( Smith et al., 2002 K 
After background subtraction had been performed the energy range was fur- 
ther reduced by discarding the energy bins with counts less than 3<r above the 
background. 
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Table 1. Flares suitable for analysis, flare positions are given in arcsecond from the Sun 
centre. 



1 


Flare Date 


Start Time (UT) 


GOES Class 


x position 


y position 




a 


20-Aug-2002 


08:25:21 


M3.4 


562 


-270 


0.72 | 


1 b 


10-Scp-2002 


14:52:47 


M2.9 


-622 


-244 


0.72 | 


c 


17-Jun-2003 


22:52:42 


M6.8 


-783 


-148 


0.52 | 


1 d 


2-Nov-2003 


17:16:00 


X8.3 


770 


-343 


0.51 | 


e 


10-Nov-2004 


02:09:40 


X2.5 


738 


116 


0.69 | 


1 f 


15-Jan-2005 


22:49:08 


X2.6 


117 


325 


0.93 | 


1 g 


17-Jan-2005 


09:43:44 


X3.8 


441 


301 


0.86 | 


1 h 


10-Scp-2005 


21:34:26 


X2.1 


-667 


-255 


0.69 | 




Figure 8. Positions of all 8 flares studied on the solar disk. The inner rings indicate 
heliocentric angles of 30° and 60° . 

4. Inversion method 

The cross-section is similarly divided into two components Q F represents the 
bremsstrahlung cross-section in the forward direction, that is the radiation beam- 
ed in the same direction as the electron was travelling and Q B represents the 
cross-section for the radiation beamed in the opposite direction to the electron. 
These matrices are determined by taking the full angular dependant cross-section 
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and averaging over a range of angles similar to the method applied to determine 
the upward and downward components of the photon flux 

Q(e,E,6 ) = 1 / Q'(e,E,P,6 )sinPdP , (12) 

2tt(1 - cos(a)) 7,3=0 

with a = 90° Q B {e, E) = Q(e, E, 6 Q = 180° -9') and Q F (e, E) = Q(e, E, 6 = 0°) 
This results in a directly observed count spectrum given by 

— \r roc 

Z o( £ ) = 4^2 J (Q F (e,E)F u (E) + Q B (e,E)F d (E))dE (13) 

and a downward directed photon flux given by 

I d(e) = £^J (Q B (e,E)F u (E) + Q F (e,E)F d (E))dE. (14) 

A Green's function approach can be used to determine the fraction of this 
downward directed photon flux reflected back towards the observer by albedo 
dKontar et al, 20060 . 

The matrix relation between the observed photon spectrum and the bi-direc- 
tional electron spectra is now given by 

C = S(Q F + AQ B Q B + A Q F )(f^)> ( 15 ) 

where A is a discretized matrix representing the Green's function and S repre- 
sents the photon to count spectral response of RHESSI. 

Equation (|T5|) must be simplified so that the array containing the cross-section 
and Green's functions is converted into a standard two dimensional matrix and 
the two component electron spectrum is represented as a one dimensional vector. 



A method of performing this transformation is described in Hubeny and Judge 
(1995])- The equation can then be solved using the standard regularised inversion 
methods. 

A significant problem is determining the electron spectrum which produced 
a given photon spectrum. In general there is a linear relationship between the 
observed count spectrum represented by vector C and the electron distribution 
which produced it, F, which can be written the form 

C = MF (16) 

here M is a matrix defined by M = SQ where Q is a matrix representing the 
bremsstrahlung cross-section in this case Q = (Q F + AQ B Q B + AQ F ) To 
determine the emitting electron spectrum this must be solved for F, here this 
is equivalent to the matrix ( ^ J . However this problem is ill-posed as noise in the 



data is amplified making a direct inversion impossible (Bertero, Mol, and Pike, 1988 



A technique for solving this is regularized inversion ( |Tikhonov, 1963] ) 
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As Equation (|16[) does not have a unique solution additional constraints 
must be imposed by considering the physical properties of the electron distribu- 
tion. The solution can then found by solving the Lagrange multipliers problem 

||MF - C|| 2 + A||LF|| 2 = minimum, (17) 

this method is known as Tikhonov regularisation. Here L represents the addi- 
tional constraint, in this case an approximation of the differential operator is used 
and A is the regularisation parameter, a variable which determines the degree of 
smoothness imposed on the solution. Equation (|17[) can be solved by Generalized 



Singular Value Decomposition for arbitrary A and L (Kontar et ai, 2004 1. 

The choice of the regularisation parameter here is determined by examining 
the normalised residuals which are defined by = ((MF)j — C^/SCk- Then 
the best value of A by minimizing 

IKmfa-cxsc)- 1 !!^"- (is) 

The choice of L is related to the constraint applied to the solution and in this 
case should be related to the underlying physics, here a finite difference matrix 
representing the first derivative is used, dependence The robustness of the solu- 



tion can be improved if the equation is first preconditioned (Kontar et al., 2004). 



A forward fit performed on dependence the data using a standard model of a 



thermal component plus a broken power law (Holman et ai, 2003). This esti- 
mated electron spectrum, Fo , is used as a starting point for the regularised 
inversion. The inversion is performed on the difference between the data, C and 
the fit, MFo, This modified data vector and the cross-section matrix are also 
scaled by a factor of yMFo- These transforms both make the solution much 
flatter and so less prone to errors. 

This method has been applied to the problem of inverting the angle averaged 
electron spectrum numerous times (e.g. IBrown et alA I2006I2005P but can be 
extended to determine an angular dependant electron flux. This is done here by 
extending the cross-section matrix to take account of the angular dependence 
of the brcmsstrahlung cross-section and then solving for an electron flux matrix 
with two components, one directed down towards the photosphere, Fj, and one 
directed towards the observer F u ( |Kontar and Brown, 2006"] ). 

The inversion algorithm is applied to the selected time intervals of the flare 
and the upward, nVF u , and downward, nVFd, electron fluxes calculated. As a 
first estimate a thermal plus double power-law fit is performed with F u = Fd- 

The errors on the electron flux components are calculated by combining the 
errors on the count flux and the errors on the background. Random perturba- 
tions are applied to the count flux based on the error and the electron flux is 
recalculated. The distribution of these realizations is then used to estimate the 
error on the electron flux. The regularised solution also has finite resolution. 
The resolution matrix is defined as R = M^ 1 M trU e, where F true = M true C is 
the true solution to the inverse problem and F\ = M^ 1 C is the regularised 
solution. The resolution matrix quantifies the horizontal errors of the solution, 
so the identity matrix (zero horizontal errors) correspond to the direct inverse 
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. For any practical situations, the regularization imposes a spread on the 
strong peaks centred on the main diagonal, this is an unavoidable occurrence in 
any inverse problem. The FWHM of each of the rows of this matrix is taken as 
the energy resolution for that energy bin and is considered here as the horizontal 
error in the electron flux (Figure 11). 

The anisotropy was defined to be the ratio of nVF^ to nVF u . Confidence 
strips for the total anisotropy were calculated using the same method as the 
errors in the electron flux (Figure 12). Random perturbations were applied within 
the "confidence river" defined by both the horizontal and vertical error bars. 

In order to test the method, we have applied it to simulated electron spectra. 
The electron spectra have been assumed to have simple functional forms and 
equivalent photon spectra and albedo reflection were added using analytical 
expressions for electron fluxes. Random noise, at a similar level to the noise 
estimated from RHESSI observations, has been added to the resulting count 
spectrum. This simulated spectrum was then inverted using the same algorithm 
as the real data. 

Spectra from isotropic initial distributions generally gives a result which is 
consistent with the input spectrum within errors and shows a reasonable dis- 
tribution of residuals. For weak anisotropy the results generally show broader 
confidence intervals than for the isotropic case with the same level of noise so 
that the solution is often consistent with the weakly anisotropic input spectrum 
and an isotropic input spectrum. For stronger levels of anisotropy (F c i/F u > 10 
at 100 keV), the method tends to give unphysical negative values for the electron 
flux. This can be avoided by increasing the regularisation parameter to force the 
solution to be smoother and ensure the solution is positive everywhere. However 
this approach leads to under-regularisation and unacceptably large residuals 
suggesting that the method cannot converge on a physically meaningful solution 
which satisfies the data. We emphasise that there is an upper limit to the size of 
the maximum anisotropy detectable as it is difficult to constrain an anisotropy 
that is greater than the fractional error in the larger component of electron flux. 

The tests also confirm that this method works best for flares with high en- 
ergy counts close to the disk centre and that the anisotropy cannot be reliably 
inferred for weak or limb events, as was expected from the forward modelling. 
All the inversions of RHESSI data show physically sound results with reasonable 
residuals. 

4.1. 20th August 2002 

This flare was detected on 20th August 2002 around 08:20 UT with the impulsive 
peak starting about 08:25 UT. It was detected with a heliocentric angle of ~ 43° 
equivalent to [i — 0.73. The flare also shows good count statistics up to 400 keV. 
As this flare had attenuator status changes from AO (open telescope) to Al (thin 
shutter in) at 08:25:16 UT and to A3 (both shutters in) at 08:25:44 UT, the anal- 
ysis was only performed over the 16 s period rather than the 64 s period studied 
for most flares. There is some particle contamination over the impulsive phase. 
This flare was extensively studied by |Kasparova, Kontar, and Brown] (I2007[) and 
the background subtraction used in this paper is similar to the subtraction 
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Figure 9. RHESSI lightcurves of each of the flares studied in 7 energy bands- black 7-12 
keV, purple 12-25 keV, blue 25-50 keV, green 50-100 kcV, yellow 100-300 keV, orange 300-800 
keV, red 800-5000 keV. The vertical lines show the accumulation time interval used. The plots 
are semi-calibrated, a diagonal approximation of the RHESSI response is used to estimate the 
photon flux from the measured counts. There are still instrumental artefacts present with the 
very sharp spikes and dips being the result of attenuator status changes. All times are in UT. 
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Figure 10. Impulsive phase count spectra accumulated by RHESSI for each flare studied 
(see Figure 9). The black line shows the background subtracted counts and the magenta line 
the background. 
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Figure 11. Results of the invertion procedure for 64 second intervals from two selected flares 
on 10th November 2004 (left) and 17th June 2003 (right). Top panel shows the measured 
count spectrum (full line) overplotted with the count spectrum corresponding to the calcu- 
lated regularised electron spectra (dashed line). The second panel shows regularised electron 
spectrum with associated l-cr vertical and horizontal error bars for each point, the light grey 
line denotes the upward electron flux and the dark grey line the downward electron flux. The 
third panel shows the normalised residuals for each time interval and the bottom panel shows 
the cumulative residuals. 



described there. This flare was previously analysed using bi-directional inversion 
bv lKontar eTaEl ((30Dg|l . 

4.2. 10th September 2002 

This flare was detected on 10th September 2002 between 14:02 and 15:15 UT 
with the impulsive peak starting about 14:52 UT. It was detected with a helio- 
centric angle of ~ 44° equivalent to \x = 0.72. The flare also shows good count 
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Figure 12. The anisotropy of the electron spectrum (defined as F^/F u ) for all eight flares 
studied. The dark grey area represents the la confidence interval and the light grey the 3cr 
confidence interval. 
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statistics up to 300 keV. This flare also has an attenuator status change from 
AO to Al at 14:52:43 UT and to A3 at 14:54:16 UT, so the impulsive phase is 
taken to be a 32 s time interval between 14:52:47 and 14:53:19 UT. 

4.3. 17th June 2003 

This flare was detected on 17th June 2003 starting at approximately 22:30 UT. 
As RHESSI shows significant particle contamination during the early stages 
of this flare analysis was performed on a later impulsive peak with accumula- 
tion starting at 22:52:42 UT. It was detected with a heliocentric angle of 59° 
equivalent to \i = 0.51. The flare also shows good count statistics up to 300 keV. 

4.4. 2nd November 2003 

This flare was detected on 2nd November 2003. It was detected with a helio- 
centric angle of ~ 59° equivalent to fi = 0.51. The flare also shows good count 
statistics up to 300 keV. RHESSI showed some elevated particle levels during 
the impulsive phase of the flare so analysis was confined to the earlier part of 
the impulsive phase. 

4.5. 10th November 2004 

This flare was detected on 10th November 2004. It was observed with a heliocen- 
tric angle of 46.5° equivalent to /j, = 0.69. The flare shows no significant particle 
measurements during the impulsive phase and low probability of pulse pileup. 
The flare also shows good count statistics up to 500 keV. 

4.6. 15th January 2005 

This flare was detected on 15th January 2005. It was detected with a heliocentric 
angle of 20° equivalent to fi = 0.93, the closest of all the flares selected to the 
disk centre and therefore the most likely to show evidence of strong downwards 
directivity. The flare also shows good count statistics up to 400 keV. 

4.7. 17th January 2005 

This flare was detected on 17th January 2005 between 09:30 and 15:15 UT with 
the impulsive peak starting about 09:42 UT. It was detected with a heliocentric 
angle of ~ 31° equivalent to \x = 0.86. The flare also shows good count statistics 
up to 300 keV. As this flare occurs in the tail of a previous flare is has very high 
background at low energies. Counts below 18 keV were not accumulated for 
this flare. This flare was also previously analysed using bi-directional inversion 
IKontar et q"D (j2"0"0"6]) . 
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4.8. 10th September 2005 

This flare was detected on 10th September 2005. It was detected with a heliocen- 
tric angle of ~ 46° equivalent to /x = 0.69. This flare showed negligible particle 
contamination and a low probability of pulse pilcup. The flare also shows good 
count statistics up to 300 kcV. Due to changes in the attenuator status during 
the impulsive phase from Al to A3 at 21:34:12 UT and the maximum time 
interval studied for this flare is 32 s starting at 21:34:26 UT. 

4.9. Short time intervals 

Four seconds is roughly the rotation period of RHESSI and thus the shortest 
time interval studied here. Longer time intervals of 8, 16, and 32 seconds were 
also studied for some flares depending on the length of the impulsive phase. 
The results for these time intervals are very similar to the results for the full 
impulsive phase. The inversions for short time intervals generally show confidence 
intervals around an anisotropy of 1 at the la level extending to around 2 below 
100 keV and sharply increasing above that. As the count statistics are lower 
for the shorter time intervals the confidence intervals are wider than for the full 
impulsive phase. There was no discernible statistically significant variation in 
the level of anisotropy for the duration of the impulsive phase. As an example 
the anisotropics for each of the 4 second time intervals for the flare on 10th 
November 2004 are shown in Figure [T21 

5. Discussion and Conclusions 

This analysis shows consistently that for almost all flares studied by RHESSI 
that the recovered Fd/F u is close to unity within the confidence intervals be- 
ing consistent with an isotropic pitch angle distribution. For almost every flare 
downward beaming of a ratio greater than ~ 3 : 1 is ruled out to 3 — a confidence 
below ~ 150 keV. The only clear exception to this is the flare on 20th August 
2002 between 30 and 50 keV, where the recovered flux appears to be inconsistent 
with isotropic at the 3-cr level and suggest a slightly greater (1.5 — 2) upward 
flux. This flare is unusual in several respects. There is a high level of particle 
contamination throughout the impulsive phase. Several background subtractions 
were examined to attempt to account for this. This flare is also one of the flattest 
flares studied, which makes pileup correction more difficult to estimate. Also, 
this is one of only two flares studied where the attenuator status was Al for the 
examined time interval. 

These measurements appear to rule out any strong beaming such as would be 
expected in the basic collisional thick target model. However as only two com- 
ponents are recovered and the confidence intervals can be fairly large using this 
method, the observations are consistent with a range of possible pitch angle dis- 
tributions including fully isotropic distributions, pancake distributions and weak 
beaming below the measured confidence level. The size of the uncertainties could 
be reduced with better count statistics and better energy resolution as each of 
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these distributions does show different spectral variation. The forward modelling 
shows very clearly that when there is significant beaming in the emitting electron 
population it will result in a very strong albedo emission. For many cases with 
substantial beaming, particularly close to the disk centre, the albedo component 
can dominate over the primary component. In order to directly compare the 
measured results to the forward model a plot of A/i against anisotropy (Fd/F u ) 
is included (Figure [T4"|) both for the functional form described by Equation 3 and 

for another commonly used form F oc exp ^ — ^jr~^J- It should be stressed that 

the results in Figures [12] and [13] are model independent and can be interpreted 
in terms of a variety of models, including, but not limited to those considered in 
Figure HU 

These results are consistent with previous published results which showed 
little evidence of directivity below 300 keV. It should be noted that this study 
measures anisotropy in terms of the electron flux, whereas for other types of 
study the parametrisation of anisotropy is often in terms of the directivity of 
the X-ray emission for stereoscopic studies and the centre-to-limb variance for 
statistical studies. These are generally related to the electron anisotropy in a 
model dependant manner. As the X-ray emission can be quite broad, particularly 
for low energies, a large anisotropy in the electron spectrum could result in a low 
photon spectrum directivity. Kasparova, Kontar, and Brown| (|2007() performed 



a centre-to-limb study using RHESSI data and inferred a directivity ratio be- 
tween 0.2 and 5 in the range 15 - 20 keV. As the emission below ~ 30 keV is 
expected to be predominantly produced by thermal electrons it is expected that 
the distribution in this energy range should be isotropic. This is particularly 
true for the flares on 17 June 2003 and 10 November 2004 which show strong 
thermal components however this may not be the case for flares which show a 
weak thermal component such as the flare on 20 August 2002. 

As this study measured X-rays in the energy range 10 - 500 keV the reli- 
ability of the inversion above approximately 250 keV is questionable. As can 
be seen from Figure [12] the confidence interval increases significantly at a few 
hundred keV. Thus it is difficult to make comparisons with the SMM studies 
which examined X-ray measurements above 300 keV. However the measure- 
ments in this study are for the most part in agreement with previous studies 
( |McTiernan and Petrosian, 1991[ ). 

As electrons propagate through the corona and chromosphere they will be 
pitch angle scattered by Coulomb collisions (e.g. ILeach and Petrosian! 119811 
|MacKinnon and Craig[ I1991[) , although it seems that collisions will be insuffi- 
cient to isotropise an initially beamed distribution ( |Brown, 1972| [Leach and Petrosian, 1981[ ) 
These results, therefore, suggest that either the accelerated electron population 
is more isotropic, or other transport effects are more important than anticipated. 
Specifically, the electron scattering by various wave-particle interactions could 
increase the pitch angle spread of the energetic electrons. Further, if the distribu- 
tion of energetic electrons is close to isotropic, the role of return current should 
be diminished. In addition we note, that although the return current itself does 
contribute to the formation of a backward going beam, it is likely to be more 
efficient at energies below ~ 50 keV, so that the higher energy electrons are 
expected to be weakly affected ( |Holman el al., 201 1[ ). 
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Figure 13. The anisotropy of the electron spectrum (defined as F^/F u ) for the first eight 4 
second time intervals for the flare occurred on 10 November 2004. The first interval starts at 
02:09:40 UT and the intervals shown here cover the most intense part of the impulsive peak. 
The dark grey area represents the lcr confidence interval and the light grey the 3<r confidence 
interval. 
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Figure 14. Pitch angle spread, Afi, for various anisotropics F^/F u using 
F(n) oc cxp ( ~ ( ^~2 )2 ) (solid line) and F(p) oc exp (dashed line). The 

vertical dotted line shows an anisotropy of 3. 
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